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Abstract 

In a recent paper ji]] we investigated Potts models on "thin" random graphs - generic Feynman di- 
agrams, using the idea that such models may be expressed as the N — ► 1 limit of a matrix model. The 
models displayed first order transitions for all q > 2, giving identical behaviour to the corresponding 
Bethe lattice. 

We use here one of the results of jlj namely a general saddle point solution for a q state Potts 
model expressed as a function of q, to investigate some peculiar features of the percolative limit q — > 1 
and compare the results with those on the Bethe lattice. 



1 Introduction and q > 2 Solutions 

In a series of papers jy, g, [| || we have exploited the observation made in Q that an ensemble of 
random graphs could be thought of as arising from the perturbative Feynman diagram expansion of a 
scalar integral to investigate spin models living on such graphs. The random graphs are of unrestricted 
topology as can be seen by thinking of the scalar model as the N — *■ 1 limit of an N x N Hermitian 
matrix model. We consequently denote them as "thin" graphs, to distinguish them from the planar "fat" 
graphs which appear in the N — > oo limit familiar from two dimensional quantum gravity. 

The random graph spin models are of interest because they display mean field behaviour due to 
the tree-like local structure of the graphs Q. Random graphs, which are closed, have an advantage 
over genuine tree-like structures such as the Bethe lattice for both numerical and analytical work because 
dominant boundary effects are absent. The numerical advantages of random graphs over the Bethe lattice 
have also been noted by other authors Q . The equilibrium behaviour on random graphs of ferromagnetic 
Ising models || and spin glasses Q was found to parallel that of the equivalent model on the appropriate 
Bethe lattice with the same number of neighbours ||. 

If we take the Hamiltonian for a q-state Potts model to be 

ff^Et^- 1 ). (!) 

<ij> 

where the spins Oi take on q values, the partition function on </> 3 (i.e. 3-regular) random graphs with 2n 
vertices is (l) 

Z»W X = i / T^T / exp(-A), (2) 
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where the action A is given by 
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K is the inverse of the quadratic form in the action, the coupling c is related to the temperature as 

1 



(exp(2/3) +q-2) 

and N n is the number of undecorated (no spin) <j> 3 graphs with 2n vertices 



(4) 
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The role of the contour integration is to pick out graphs with 2n vertices and the prefactor of N n 
disentangles the factorial growth in the number of graphs themselves from any non-analyticity due to a 
transition in the decorating spin model. 

The model may be solved in a saddle point approximation, valid for large n, and one finds a high 
temperature solution of the form <pi — 1 — (q — l)c, Vi and low temperature broken symmetry solutions 
<pi — . . . 4> q -i 7^ 4> q . There are two low temperature branches, the first being 



1 - (g - 3)c - y/1 - 2(g - l)c+ (q - 5)(g - 1)1 

h...q-l - r 



1 + (g - l)c + y/1 - 2(q - l)c + (q - 5)(g - l)c 2 

the second having the signs in front of the square roots reversed 

1 - (<? - 3)C + ^l-2(q-l)c+(q-5)(q-l)c2 
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1 + ( 9 - l)c - VI - % - l)c + (9 - 5)(g - l)c 2 
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The low temperature solutions can be recovered from the saddle point equations for an "effective 
action" A, namely dA/d<fi = dA/ddj) = where 

A = \{q - 1) [1 - c(q - 2)] fa 1 - l(q - l)^ 3 + ^ 2 - ^ 3 - c(q - 1)# (8) 

with <fi = ^i... g _i, <p — 4> q . Similarly, at high temperature one can write another effective action 

4) = f(l-c(g-l))$-|$ (9) 

whose saddle point equation gives fa = 1 — (q — l)c. Both equs.(||,||) follow from imposing the expected 
symmetry on the fa in the original action in equ.(^). 
The Potts magnetisation 

m = mbm (10) 

is related to the standard order parameter by 

M = gmax(m) ~ 1 (11) 

q - 1 

and gives a clear picture of the topology of the phase diagram. In Fig.l we plot m against c for q = 4, 
this being representative of all the q > 2 models. The upper (solid) branch is given by m for the solutions 
in equ.(^|) where it gives a maximum and the lower (dashed) branch is given by m for the solutions in 
equ.(0). The high temperature (dotted) line is at m = m = 1/q. 
A first order transition occurs for q > 2 at 

c(Q) = ■ (12) 

q z 

when the low temperature and high temperature saddle point effective actions become equal. One finds 
a jump in the magnetisation at the first order transition point of 

AM = ^— -. (13) 
q - 1 

Spinodal points are present at P where the high temperature solution joins the lower branch 

and O where the low temperature branches first become real 

C(0) = (*-D(* -5) ■ (15) 

For q = 2 0,P,Q merge and one recovers the continuous mean field Ising transition. All the exponents, 
critical values and jumps are in agreement with the calculations in M] for the Bethe lattice proper. 



2 Percolation and the q — ► 1 limit 

As we have a solution and expressions for the critical points that are apparently valid for all q we shall 
now attempt to emulate the work of (ll[] for the Bethe lattice and explore the percolative limit q—*l. The 
solution of the Bethe lattice percolation problem is of course well known from elementary considerations 
Jl2] , ^3). One finds that the percolation threshold on a Bethe lattice with co-ordination number z is 
p c = l/(z — 1) and other cluster quantities are also explicitly calculable. Our objective here, as in |[l| 
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for the Bethe lattice proper, is rather to shed light on the properties of the Potts model correspondence 
with percolation. 

The mapping between the Potts model and a percolation problem was first written down by Fortuin 
and Kastelyn |H and is the basis of the cluster algorithms that have been so successful in combating 
critical slowing down in simulations ]lq| . The q-state Potts model partition function may be written as 
a bond partition function in terms of spin clusters 

z = E p°(i-p) u q Nc (i6) 

con figs 

where o = fl occupied bonds, u = jj unoccupied bonds and N c = § clusters. The factor p, which converts 
the "geometrical" spin clusters into critical spin clusters is 1 — exp(— 2(3) and the sum over configurations, 
in distinction to standard lattices, includes a sum over different Feynman diagrams. From equ.(|l6|) 
one can see that the mean number of connected clusters per vertex N a (averaged over the ensemble of 
Feynman diagrams) is given by 

No = Um rim -^S(Z n ) 

n^oo q~*l n(q — 1) 

where Z n is the partition function in equ.(|^). Other thermodynamic quantities in the Potts model are 
related to percolative quantities, for instance the magnetisation m to the percolation probability P 

and the magnetic susceptibility \ to the mean size of finite clusters S 

S (1 - P) = Hm —X (19) 
g-i (q - 1) 

In principle the above relations apply per site, but with the uniform random graphs we consider it is 
possible to consider spatially averaged values and drop any site indices. 

To get some understanding of the limit q — » 1 that is required in the above let us first look at the 
region 1 < q < 2. In Fig. 2 we have plotted the values of c(0),c(P) and c(Q) where we can see that 
they fan out again from equality at q = 2 in the same order as for q > 2 (i.e. Q is sandwiched between 
an upper value at O and a lower value at P). We can also see that c(0) and c(Q) diverge as q — > 1. 
The sign of the jump in the magnetisation is reversed with respect to the q > 2 solutions as can be seen 
in Fig. 3, where we have taken q = 1.2 as an illustrative example. On reducing c (i.e. T) a first order 
transition occurs at Q from the horizontal high temperature solution to the lower dashed branch of the 
low temperature curve. The portions QP and the upper branch to the left of P are metastable, whereas 
the portion OQ of the lower branch and the remainder of the high temperature line to the left of P are 
unstable. 

As we have noted, as q — > 1, c(Q) — > oo or C rit ~~ > 0- The first order transition is therefore clearly 
not that associated with the percolation problem as 

p = 1 - exp(-2/3) (20) 

and the known Bethe lattice solution (for z = 3) is p cr it — 1/2. We must look rather at the spinodal 
point P where one has c(P) — l/(2q — 1), or exp(2/3(P)) = q + 1. In the limit q — > 1, one thus has 
exp(2/3(P)) — > 2 and p CT it —* 1/2, as expected. This result is not a surprise, given the discussion of [0 
for percolation on the Bethe lattice, where an exactly analogous situation occurs. In the limit q — > 1 the 
percolation transition is thus associated not with the first order transition at Q to a stable state, but 
rather to the spinodal point P where the high temperature phase joins a metastable branch. 

It is also possible to directly calculate various percolative quantities in the saddle point approximation. 
For example, to find the percolation probability P we note <f> — > c, <f> — > 1 as q — > 1, so we have from 
cqus. 

(|0j|li) P = 1 -c 3 , which using the relation between c and p for q = 1, namely c = (1 — p)/p, gives 



/' 1 ( — ) • (21) 
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This reproduces the standard result for P on a Bethe lattice with 3 neighbours. The probability of a 
point being connected to infinity by occupied sites is pP with our conventions. 

More difficult is the calculation of the susceptibility, which requires the solution of the saddle point 
equations for the effective action in an external field H 



A= l -{q-l)[l-c{ q -2)]4> 2 -\(q-l) 



I 7 



H ~ 
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(note the H in front of <fi 3 ), as we have 
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dH Y dH 



+ ^3 



(22) 



(23) 



The first two terms in equ.(|23|) come from the implicit dependence of the solutions on H and the final 
term is dm/dH. In calculating this one uses the expression for the magnetisation in the presence of a 
field 

m = ^ (24) 

{q-l)^ +H^ 

The full solutions for 4>, cj> in the presence of a field are not particularly illuminating so wc do not reproduce 
them here, but we find the simple limits 

lim lim d> = c 

ff->l 



lim lim — — 

OH 

lim lim d> 

lim lim — — 

q ->l h^i dH 



giving 



lim-^- = Sil-P) 
a-*i q—1 

which again reproduces the Bethe lattice result for S. 



c-1 
1 



= -1 



;1+2C 



(25) 



(26) 



3 Discussion 

The somewhat peculiar features of the correspondence between the q — * 1 limit of the Potts model and 
percolation on the Bethe lattice in which the spinodal point is tied to the percolative transition are 
clearly preserved on Feynman diagrams. Using the relations between Potts model observables and those 
associated with the percolation problem, we obtain the standard Bethe lattice results for quantities such 
as P and S. The inventory of results in which the (ferromagnetic) transitions of spin models on random 
graphs are identical to those on corresponding Bethe lattices is thus extended to percolation and Potts 
models analytically continued to real values of q. 

It is also interesting to note that the Ising spin glass transition on random graphs, which appears in 
the k — > limit of a k Ising replica model, displays some very similar features to the percolative limit 
discussed here. In such models the spin glass transition in the (quenched) k — * limit is associated 
with a continuous transition which appears V7c > 2 between the high temperature phase and a phase 
that is metastable when k > 2 ||. The "true" transition in the k Ising replica model is first order for 
k > 2 |u| H and is not that associated with the k — * limit. One might hope to be able to analyse 
the multi-Ising models in a similar fashion to the Potts models by writing down an effective action for 
various k and seeing if a general formula could be extracted to allow analytical continuation in k. The 
superficial similarity with the percolative limit is intriguing and certainly merits further investigation to 
see if there is any deeper connection. 
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Figure 1: The magnetisation m for a 4 state Potts model as calculated from the saddle point solutions. 
The high temperature branch is shown dotted, the upper low temperature branch solid and the lower low 
temperature branch dashed. 
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Figure 2: c(Q) (solid middle line), c(0) (upper dotted line) and c{P) (lower dashed line) plotted against 
q for 1 < q < 2. c(0) and c(Q) diverge as q — > 1. 
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Figure 3: The magnetisation m for a q — 1.2 state Potts model as calculated from the saddle point 
solutions. The key is as for Fig.l, but note that the first order transition at Q is now to the lower branch. 
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